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The supersymmetric CP problem is studied within superstring-motivated exten- 
sions of the MSSM with an additional U(l)' gauge symmetry broken at the TeV scale. 
This class of models offers an attractive solution to the /i problem of the MSSM, in 
,^Zh ■ which U(l)' gauge invariance forbids the bare jx term, but an effective \i parameter 

is generated by the vacuum expectation value of a Standard Model singlet S which 
has superpotential coupling of the form SH u Hd to the electroweak Higgs doublets. 
The effective [i parameter is thus dynamically determined as a function of the soft 
supersymmetry breaking parameters, and can be complex if the soft parameters have 
nontrivial CP-violating phases. We examine the phenomenological constraints on the 
reparameterization invariant phase combinations within this framework, and find that 
the supersymmetric CP problem can be greatly alleviated in models in which the 
phase of the SU(2) gaugino mass parameter is aligned with the soft trilinear scalar 
mass parameter associated with the SH u Hd coupling. We also study how the phases 
filter into the Higgs sector, and find that while the Higgs sector conserves CP at the 
renormalizable level to all orders of perturbation theory, CP violation can enter at the 
nonrenormalizable level at one-loop order. In the majority of the parameter space, 
the lightest Higgs boson remains essentially CP even but the heavier Higgs bosons can 
exhibit large CP-violating mixings, similar to the CP-violating MSSM with large \x 
parameter. 



1 Introduction 



While phenomenological models with low energy supersymmetry (SUSY) are arguably the 
best candidates for physics beyond the Standard Model (SM), they typically include a large 
number of parameters associated with the soft supersymmetry breaking sector. For example, 
the minimal supersymmetric standard model (MSSM), which has two Higgs doublets and 
conserved R-parity, contains 105 new parameters [1], including the bilinear Higgs superpo- 
tential parameter \i and the soft SUSY breaking parameters (this counting does not include 
the gravitino mass and coupling). The parameter count generically increases if such SUSY 
models are extended beyond the minimal gauge structure and particle content of the MSSM 
(unless symmetry relations exist in the theory which relate subsets of parameters). Many 
of these new parameters are phases, which both provide new sources of CP violation and 
modify the amplitudes for CP-conserving processes. Even if certain sectors of the theory 
exhibit no CP violation at tree level (e.g. if the relevant phases can be eliminated by global 
phase rotations), the phases can leak into such sectors of the theory at the loop level and 
have an impact on collider phenomenology and cosmology. 

In contrast to the SM, in which the only source of CP violation is present in the CKM 
matrix and thus is intimately tied to flavor physics, CP-violating phases within SUSY mod- 
els can occur in both flavor-conserving and flavor- changing couplings. The phases of the 
flavor-conserving couplings (which have no analogue in the SM) are of particular interest be- 
cause they can have significant phenomenological implications which can be studied without 
knowledge of the origin of intergenerational mixing. In the MSSM, these phases are given 
by reparameterization invariant combinations of the phases of the gaugino mass parameters 
M a (a = 1, 2, 3), the trilinear couplings Af, the // parameter, and the Higgs bilinear coupling 
b = fiB. However, not all of these phases are physical; after utilizing the U(1)pq and U(1)r 
global symmetries of the MSSM, the reparameterization invariant phase combinations are 
@f = + ^A-s ~~ and 9 a = <f)^ + <pM a ~ <Pb (in self-evident notation). 

Such phases have traditionally been assumed to be small due to what is known as the 
supersymmetric CP problem: the experimental upper limits on the electric dipole moments 
(EDMs) of the electron, neutron, and certain atoms individually constrain the phases to 
be less than 0(1CT 2 ) (for sparticle masses consistent with naturalness) [2, 3, 4]. However, 
recent studies have shown that EDM bounds can be satisfied without requiring all repa- 
rameterization invariant phase combinations to be small, if either (i) certain cancellations 
exist between different EDM contributions [5, 6, 7, 8, 9], or (iii) the sparticles of the first 
and second families have multi-TeV masses [10]. 1 Even within each of these scenarios, the 
particularly strong constraints arising from the atomic EDMs [13] lead to a general upper 
bound of < O(10~ 3 ) on the reparameterization invariant phase present in the chargino sec- 
tor (#2 in our notation), while the other phases are comparatively unconstrained [8]. These 
constraints will be discussed in detail later in the paper; for now, it is worth noting that 

1 Thc EDM bounds are more difficult to satisfy in both of these scenarios when tan/3 (the ratio of 
electroweak Higgs VEVs) is large. Not only are cancellations in the one-loop EDMs more difficult to achieve, 
but certain two- loop contributions are then enhanced [11, 12] which do not decouple when the sfcrmions are 
heavy. In part for this reason, we will restrict our attention in this paper to the small tan/3 regime. 
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this "CP hierarchy problem" is an intriguing issue to be addressed within models of the soft 
parameters which include CP violation. 2 

Of course, if the reparameterization invariant phases are sizeable, they can have important 
phenomenological consequences. Within the MSSM, one of the examples in which these 
phases can have a significant impact is the Higgs sector. As is well known, the MSSM Higgs 
sector conserves CP at tree level. However, radiative corrections involving the SM fields and 
their superpartners (with the dominant effects typically due to top and stop loops) have 
a substantial impact on Higgs masses and mixings. For example, the one-loop radiative 
corrections substantially elevate the tree-level theoretical upper bound of Mz on the mass 
of the lightest Higgs boson [16]; these results have been improved by utilizing complete 
one-loop on-shell renormalization[17], renormalization group methods [18], diagrammatic 
methods with leading order QCD corrections [19], two-loop on-shell renormalization [20], 
and complete two-loop effective potential [21]. Indeed, if the radiative corrections include a 
nontrivial dependence on phases, the Higgs potential violates CP explicitly at one-loop. The 
Higgs mass eigenstates then no longer have definite CP properties, which leads to important 
implications for Higgs production and decay [22, 23, 24, 25]. 

The MSSM offers a minimal framework for stabilizing the Higgs sector against quadratic 
divergences. However, it is well known that the MSSM has a hierarchy problem with respect 
to the scale of the superpotential fx parameter [26], which has a natural scale of O(Mgut), 
and the electroweak scale. An elegant framework in which to address this 11 fx problem" is to 
generate the fx parameter via the the vacuum expectation value (VEV) of a SM singlet S. 
One simple possibility 3 [28] is to invoke an additional nonanomalous U(l)' gauge symmetry 
broken at the TeV scale, as expected in many string models. For suitable U(l)' charges, the 
bare fx parameter is forbidden but the operator h s SH u -H d is allowed, such that an effective fx 
term is generated after S develops a VEV of order the electroweak/TeV scale (assuming the 
Yukawa coupling h s ~ 0(1), as is well-motivated within semirealistic superstring models). 
This framework is of particular interest because such extra U(l) groups are often present in 
plausible extensions of the MSSM, and in fact are ubiquitous within many classes of four- 
dimensional superstring models. Indeed, additional nonanomalous U(l) gauge groups are 
present in virtually all known 4D string models with semirealistic features, such as gauge 
structure which includes SU (3) c x SU (2) L x U(1) Y (or a viable GUT extension) and particle 
content which includes the MSSM fields. 4 

2 However, there are unavoidable theoretical uncertainties involved in the determination of the hadronic 
EDMs and the atomic EDMs (see e.g. [14, 15] for discussions). These uncertainties are particularly problem- 
atic for the mercury EDM, which yields the strongest constraints on the SUSY phases. For this reason, there 
are disagreements in the literature over how to include this bound. Here we take a conservative approach 
by including the Hg EDM constraint. 

3 Thc fi parameter can also be generated in models with no additional gauge groups, i.e. the next-to- 
minimal supersymmetric standard model (NMSSM). However, NMSSM models generically possess discrete 
vacua and the tensions of the walls separating them are too large to be cosmologically admissable [27]. 

4 For example, many examples of such semirealistic models have been constructed within perturbative 
heterotic string theory (see e.g. [29] for an overview). An interesting class of constructions is the set of 
free fermionic models [30, 31, 32], in which a number of extra Z7(l)'s are always present at the string scale. 
Whether or not all of these U(l)s persist to the TeV scale depends on the details of the vacuum restabilization 
procedure. Although there are cases in which only the MSSM gauge structure remains at low energy [33], 
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Within this class of models, the electroweak and U(l)' symmetry breaking is driven by 
the soft SUSY breaking parameters, and hence the Z' mass is expected to be of order a few 
TeV or less. Such a Z' should be easily observable at either present or forthcoming colliders. 
Indeed, the nonobservation to date of a Z' puts interesting but stringent limits on the Z' 
mass and mixing with the ordinary Z both from direct searches at the Tevatron [38] and 
indirect tests from precision electroweak measurements [39]. Although limits depend on the 
details of the Z' couplings, typically M%i > 500 — 800 GeV and the Z — Z' mixing angle 
a z~z r O(10~ 3 ). 5 These models have been analyzed at tree level in [41, 42, 43], where it 
was found that there are corners of parameter space in which an acceptable Z — Z' hierarchy 
can be achieved. Further studies of a different class of string-motivated U(l)' models can be 
found in [44]. 

As the phase of the \x parameter filters into the amplitudes for many physical observables 
in the MSSM (and plays an important role in the Higgs sector at one-loop), it is worthwhile 
to analyze models which solve the fi problem in the presence of explicit CP violation. In 
this paper, we thus study the supersymmetric CP problem in U(l)' models, focusing on the 
radiative corrections to the Higgs sector of the U(l)' model of [41] in the case that the soft 
supersymmetry breaking parameters have general CP-violating phases (radiative corrections 
in the CP-conserving case has been studied in [45]). We begin by classifying the reparame- 
terization invariant phase combinations and comment on the phenomenological constraints 
on these phases from EDM bounds. We then turn to the Higgs sector, which conserves CP 
at tree level, but phases enter the Higgs potential through the stop mass-squared matrix 
at one-loop (just as in the MSSM). The VEV's of the electroweak Higgs doublets H u ^ and 
singlet S are then determined by minimizing the loop-corrected Higgs potential. Within 
this framework, an effective \x parameter of the correct magnitude is generated which also 
has a phase governed by the phases of the soft SUSY breaking parameters. We study the 
pattern of Higgs masses and mixings including the EDM and Z' constraints, and discuss the 
phenomenological implications for Higgs searches. 

2 The SUSY CP Problem in U{\)' Models 

We study the class of U(l)' models of [41], in which the gauge group is extended to 

G = SU(3) C x SU(2) L x U(l) y x U(l)', (1) 

with gauge couplings g 3 , g 2 , gy, gv, respectively. The matter content includes the MSSM 
superfields and a SM singlet S, which are all generically assumed to be charged under the 

typically one or more extra U(l)s persists to the electroweak scale [34, 35]. Additional U (l)s also are generic 
in supersymmetric braneworld models derived from Type II string orientifolds [36] (due at least in part to 
the U(N) gauge groups associated with a stacks of D branes). Phenomenological analyses also indicate that 
typically extra ?7(l)s are present in the low energy theory and broken at the electroweak/TeV scale [37]. 

5 A potentially more stringent limit on the Z' mass arises from cosmology if the U(l)' gauge symmetry 
forbids the standard implementation of the seesaw mechanism for neutrino masses. In such scenarios, the 
right-handed neutrinos may be light, and BBN constraints then require modcl-dcpcndcnt limits that in some 
cases are as strong as Mz> 4 TeV [40]. 
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additional U(l)' gauge symmetry. Explicitly, the particle content is: ~ (1, 2, — 1/2, Ql), 
E? ~ (l,l,l,g s ), 4 ~ (3,2,1/6,Q Q ), t/F ~ (3,1, -2/3, Qt/), D\ ~ (3, 1, 1/3, Q D ), if d ~ 
(1,2, -1/2, Q d ), 7? u ~ (1,2, 1/2, Q tt ), 5 ~ (1,1,0,Q 5 ), in which % is the family index. 6 

The superpotential includes a Yukawa coupling of the two electroweak Higgs doublets 
H ud to the singlet S, as well as a top quark Yukawa coupling: 

W = h s SHu ■ H d + h t U'Q 3 ■ H u . (2) 

Gauge invariance of W under {7(1)' requires that Q u + Qd + Qs — and Qq 3 + Qu 3 + Q u = 0. 
This choice of charges not only forbids the "bare" \x parameter but also a Kahler potential 
coupling of the form H u H d + h.c. required for the Giudice-Masiero mechanism [47] (the 
Kahler potential is otherwise assumed to be of canonical form). 7 

The form of (2) is motivated by string models in which a given Higgs doublet only has 0(1) 
Yukawa couplings to a single (third) family. We will consider the small (H u ) / (H d ) = tan (3 
regime only 8 such that the Yukawa couplings of the b and r can be safely neglected. The 
origin of the Yukawa couplings of the first and second generations of quarks and leptons 
is not addressed. As we are primarily interested in the third family, we shall suppress the 
family index in what follows. 

The soft supersymmetry breaking parameters include gaugino masses M a (a — 1, 1', 2, 3), 
trilinear couplings A s and A t , and soft mass-squared parameters m 2 : 

~£soft = (E M aKK + A s h s SH u ■ H d + A t h t U c Q ■ H u + h.c.) + m 2 u \H U \ 2 + m 2 d \H d \ 2 

a 

+ m^| 2 + M||Q| 2 + Mg|^| 2 + Mj|i5| 2 + M||E| 2 + M||Z| 2 . (3) 

These soft SUSY breaking parameters are generically nonuniversal at low energies. We do 
not address the origin of these low energy parameters via RG evolution from high energy 
boundary conditions in this paper. 

The gaugino masses M a and soft trilinear couplings A s>t of (3) can be complex; if so, they 
can provide sources of CP violation (without loss of generality, the Yukawa couplings h s>t can 
be assumed to be real). However, not all of these phases are physical, just as the case in the 
MSSM. Let us first consider the MSSM. The reparameterization invariant combinations of 
phases in the MSSM are easily determined by forming invariants with respect to the global 
U(1)pq and U(1)r symmetries present in the limit that the soft breaking parameters and 
the \x term are set to zero [49]; for reference, the U(1) R:P q charge assignments are presented 

6 Note that if the f7(l)' charges are family nonuniversal they provide a tree-level source of FCNC. Phc- 
nomenological bounds thus dictate that the charges of the first and second families should be identical to 
avoid overproduction of FCNC without fine-tuning [46] . 

7 Other than these constraints, we prefer to leave the U(l)' charges unspecified because our aim is not 
to construct a specific model. In an explicit model there will be additional constraints on the U(l)' charges 
(e.g from anomaly cancellation). We also do not consider kinetic mixing in the analysis [48]. However, even 
if kinetic mixing is absent at tree level will be generated through 1-loop RG running if TrQyQv ^ [48]. 

8 Here low values of tan/3 such as tan/3 = 1 are allowed (this region is excluded in the MSSM). The reason 
is that the Higgs bosons are generically heavier in U(l)' models (as in the NMSSM and other models with 
extended Higgs sectors), and even at tree level the lightest Higgs boson can easily escape LEP bounds. 
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Table 1: The U(1)r, pq charge assignments for the MSSM fields and spurions. 
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Table 2: The U(1)a charge assignments for the fields and spurions in the U(l)' framework. Note that S 
(whose VEV induces an effective \i parameter) has nonzero R charge. 

in Table 2. A convenient basis of the resulting reparameterization invariant phases thus is 
9f — ( t ) ^ J r ( t ) A f — <i ) b and 9 a = (p^+^Ma —<Pb, which enter the mass matrices of the sfermions and 
the gauginos/Higgsinos, respectively. An analysis of the MSSM tree level Higgs sector also 
suggests it is useful to exploit U(1)pq to set 0^ = (06 is then dropped from the invariants 
above), in which case the Higgs VEVs are real. 

Performing the same exercise in the U(l)' framework, one immediately notices that the 
U(1)pq symmetry of the MSSM is embedded within the U(l)' gauge symmetry. However, a 
nontrivial U(l)p symmetry remains; the U(1)r charges of the superfields and the associated 
spurion charges of the soft parameters are presented in Table 2. The reparameterization 
invariant phase combinations are therefore 9ff = 4>A f — (f>A fl , f , 9a f — <pM a — <f>A f , and 
9 ah = <pM a — 0m 6 j o , of which only two are linearly independent (e.g. <p a b = <f> a f — <fib /)• We 
will see that (in analogy to the MSSM) the tree-level Higgs sector suggests it is convenient to 
measure all phases with respect to the phase of A s . (In fact, one can go further and exploit 
the U(1)r symmetry to set (pA s = 0, although we prefer not to do that in this paper.) A 
basis of reparameterization invariant phase combinations can then be chosen as 

Ofs = <Pa } -<Pa s 

9 a s = <pM a -<pA 3 - (4) 

To see this more explicitly and to lay the foundation for our analysis of the Higgs sector 
including one-loop radiative corrections, let us now review the tree-level Higgs potential 
analyzed in [41] . Gauge symmetry breaking is driven by the VEVs of the electroweak Higgs 
doublets H u , H d 

and the singlet S. The tree level Higgs potential is a sum of F terms, D terms, and soft 
supersymmetry breaking terms: 

Vtree = V F + V D + V soft , (6) 
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in which 



V F 



\H U ■ H d 



+ \S\ 2 (\H u \ 2 + \H d 



(7) 



Vd = ^(\H u \ 2 -\H d \ 2 ) 2 + f(\H u \ 2 \H d \ 2 -\H u .H d \ 2 ) 
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+ 9 -f(Qu\H u \ 2 + Q d \H d \ 2 + Q s \S\ 2 ) , (8) 

= m 2 u \H u \ 2 + m 2 \H d \ 2 + m 2 s \S\ 2 + (A s h s SH u ■ H d + h.c). (9) 

where G 2 = g\ + g\ and gy = ^3/5gi, g\ is the GUT normalized hypercharge coupling. 
At the minimum of the potential, the Higgs fields are expanded 9 as follows: 

(S) = -^e^ivs + ^ + t^), (10) 

in which v 2 = v 2 u + v d = (246 GeV) 2 . In the above, a phase shift e td has been attached to 
(S). Since gauge invariance dictates that only the phase of the combination SH U ■ H d enters 
the potential, we can assume that the VEVs of H u d are real and attach a phase only to S 
without loss of generality (this choice is also consistent with our assignment of U{1)r charges 
in Table 2). The effective \i parameter is generated by the singlet VEV (S): 

h s v s i$ /i 1 \ 

Ms = e • (n) 

The only complex parameter which enters the Higgs potential at tree level is A s . However, 
the global phases of the Higgs fields (more precisely, of the combination SH U ■ H d ) can always 
be chosen to absorb the phase <pA s of A s by performing a U(1)r rotation on the fields, such 
that A s and the VEVs can all be taken to be real without loss of generality [41]. To state 
this another way, the minimization conditions with respect to the CP odd directions v?i,2,s 
all lead to the condition 

sin(0 + A J = O, (12) 

such that 9 = —4>a s at tree level. With this condition, the Higgs sector is CP conserving. 
The Higgs mass eigenstates thus have definite CP quantum numbers, with three CP even 
Higgs bosons -£^=1,2,3 and one CP odd Higgs boson A , as well as a charged Higgs pair 
H ± . Expressions for their masses at tree level and a discussion of the associated Higgs 
phenomenology can be found in [41]. 

9 As discussed in [41], gauge rotations can be used to set = 0. However, (HJ) = is not automatic 

and imposes constraints on the parameter space of the model. Indeed, (HJ) — implies that the physical 
charged Higgs is nontachyonic (M^ ± > 0). 
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Although the Higgs sector conserves CP at tree level whether or not the soft SUSY 
breaking parameters are complex, this is generically not true for other sectors of the theory 
and care must be taken in the phenomenological analysis (e.g. for the EDM bounds) if there 
are nontrivial CP- violating phases in the soft terms even if the Higgs sector is only analyzed 
at tree level. Clearly, this is due to the fact that the phases which enter the mass matrices of 
the sfermion and the gaugino/higgsino sectors involve the phase of the singlet VEV 9 (i.e., 
the phase of the effective \i parameter \i e fi) as well as the phases of the A terms and the 
gaugino masses. For example, the reparameterization invariant phase combination which 
enters the chargino mass matrix within this class of U(l)' models is 

®X ±= + = 4>M 2 - 4>A S + ■ ■ ■ = 02s + ■ ■ ■ , (13) 

in which the terms represented by (+ . . .) are higher-loop contributions. As previously men- 
tioned, this phase is strongly constrained by EDM experimental bounds (although the precise 
constraints can depend in detail on the other parameters of the model). More generally, the 
statement of (flavor-independent) SUSY CP violation within U(\)' models is that if any 
of the phases 9f s and 9 as defined in (4) are nonzero, they can lead to CP- violating effects 
which may be in conflict with experiment and must be checked. This is in direct analogy 
to the statement of flavor-independent SUSY CP violation in the MSSM. However, as /i is 
dynamically generated within U(l)' models, its phase M is now a function of the phases of 
the other soft breaking parameters rather than an independent quantity. 

Returning now to the question of CP violation in the Higgs sector, (12) suggests that it 
is natural to consider the combination of phases 

9 = 9 + <f> As (14) 

as the parameter which governs CP violation in the Higgs sector. Note that 9 is a reparame- 
terization invariant quantity, while 9 is not (9 = 9 in the basis in which = 0). While 9 = 
at tree level, it acquires a nonzero value at one-loop if the sfermion and gaugino/higgsino 
mass matrices have nontrivial phases. This calculation is outlined in the next section. 

3 Higgs Sector CP Violation 

Previously we discussed the SUSY CP problem within U(l)' models, and reviewed the tree 
level Higgs sector (the patterns of gauge symmetry breaking which led to an acceptable 
Z — Z' hierarchy were analyzed at tree level in [41]). In what follows, we will compute the 
one-loop radiative corrections to the Higgs sector of this class of Z' models within a general 
framework including nontrivial CP violation (radiative corrections in the CP-conserving case 
were previously presented in [45]). 

3.1 Radiative Corrections to the Higgs Potential 

The effective potential approach provides an elegant way of determining the true vacuum 
state of a spontaneously broken gauge theory. The potential has the form 

V = V tree + AV + . . . , (15) 
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where V tree is defined in (6), and the one-loop contribution AV has the Coleman- Weinberg 
form 

1 L . „ M 2 (H u ,H d ,S)_ _ 3\j 



64tt 2 



StrM 4 {H u ,H d ,S) In 



Q 2 



(16) 



in the mass-independent renormalization scheme DR. 10 In the above, Str = J2j{~ 1) 2J+1 (2J+ 
1) is the usual supertrace, Q is the renormalization scale, and M represents the Higgs field- 
dependent mass matrices of the particles and sparticles of the theory. 11 

Here we will include only the dominant terms due to top and scalar top quark loops: 



Ay 



6 



64tt 2 



£ K) ! 



k=l,2 



mA 



In 



Q 2 



-2(m t 2 ) S 



In 



V 2 



(17) 



in which the masses depend explicitly on the Higgs field components (note that Ay naturally 
vanishes in the limit of exact SUSY). The top mass-squared is given by m 2 = h 2 \H u \ 2 , and 
the stop masses-squared are obtained by diagonalizing the mass-squared matrix 



M 2 = 



M 2 LL 

M LR 



M 2 LR 
M 2 RR 



(18) 



via the unitary matrix S t as S\M 2 S t = diag (ra~ ,m~). The entries of M 2 are given by 



M 



LL 



M RR 



M LR 



= M^ + h 2 \H, 




-f) (\H U \ 2 - \H d \ 2 ) 
+ 9y'Qq [Qu\H u \ 2 + Q d \H d \ 2 + Q S \S\ 2 ^ 
= M 2 ~ + h 2 t \H u \ 2 - l -g 2 Y {\H u \ 2 -\H d \ 2 ) 

+ 9y'Qu (Qu\H u \ 2 + Q d \H d \ 2 + Q S \S\ 2 ^ 
= h t (A* t H° u *-h s SH° d ), 



(19) 

(20) 
(21) 



in which we have emphasized the fact that the LR entry depends only on the neutral compo- 
nents of the Higgs fields. As the stop LR mixing can be complex, the term hth s A t SH®H d + 
h.c. present in (M^l can provide a source of CP violation in the Higgs sector. From the 
discussion of the previous section, we can infer that this source is the phase 9 ts = <f)A t ~ 'pAs- 



10 See Martin's paper in [21] for a detailed discussion of the regularization and renormalization scheme 
dependence of the effective potential. 

11 While the complete effective potential is scale invariant, it is scale dependent when truncated to any finite 
loop order in perturbation theory due to the renormalization of the parameters and the Higgs wavefunctions. 
In the MSSM, most of the scale-dependent terms can be collected in the pseudoscalar mass, which itself can 
be regarded as a free parameter of the theory. The remaining Q 2 -dependence arises from the D term 
contributions generated by wavefunction renormalization, such that in the limit in which g-i = gy = 
all of the scale dependence can be absorbed into the pseudoscalar mass. For the U(l)' models, the scale 
dependence can be absorbed into the pseudoscalar mass only if the D term contributions vanish and the 
superpotential parameter h s — 0, because the potential also includes quartic Higgs couplings which arise 
from F terms. These properties are manifest in the Higgs mass-squared matrix presented below. 
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The vacuum state is characterized by the vanishing of all tadpoles and positivity of the 
resulting Higgs boson masses. Recalling the expressions for the Higgs fields in (10), the 
vanishing of tadpoles for V along the CP even directions <f) U) d, s enables the soft masses m^, 
to be expressed in terms of the other parameters of the potential: 





= Ml 


cos 2 /? 


-Kv 2 u 


ml 


= Ml 


sin 2 /3 


- ^vj 


m 2 s 


= Ml 


cot 2 a 





1 fdAV\ 



Vd V d( t>d /, 
1 ( dAV \ 

v s \ d(p s ) y 



(23) 
(24) 



o 

in which the subscript indicates that the derivatives of AV are to be evaluated at 0« = 
and ^ — 0. Here the various A coefficients represent the quartic couplings in the potential 

11 1 2 

K,d = -G + -Q U)d g Y > , = -jQ s 9y> > 

^ud = ~~£* 2 + QuQd 9y> + h 2 s , A USj(is = Q s Qu,d 9y' + • (25) 

The Higgs soft masses (22) are written in terms of two angle parameters: (i) tan/5, which 
measures the hierarchy of the Higgs doublet VEVs, and (ii) cot a = (i>sin/3cos/3)/i> s , which 
is an indication of the splitting between the U(l)' and electroweak breaking scales. For 
convenience, we have also introduced the mass parameter 

„ , 9 h s \A s \v s cos9 . . 

M ° = — a' (26) 

V 2 sin p cos p 

which corresponds to the mass parameter of the CP odd Higgs boson of the MSSM after 
using the definition of the effective \i parameter in (11). 

While the vanishing of the tadpoles along the CP odd directions ip Utd , s led to (12) at tree 
level, once the loop corrections are included they lead to the following conditions: 

M2sin/3cos/3tan0 = (27) 

v d V d(p u J 

A/2 sin /3 cos /3 tan = — (^—) (28) 

v u V d( Pd J o 

M 2 cotatan0 = , (29) 

v s V dip s J 

demonstrating explicitly that the phase 9 = 9 + 4>a s associated with the phase 9 of the singlet 
VEV is indeed a radiatively induced quantity. Indeed, the derivatives of AV^ with respect to 
<Pu,d,s are nonvanishing provided that there is a nontrivial phase difference between A t and 
A s (i.e., if 9 ts ^ 0). In fact, (27)-(29) all lead to the same relation for 9: 

sm9 = -f3 ht rj\sm9 ts T{Q 2 ,m\,m\) (30) 

I -"-s I 
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in which (3 ht = 3h 2 /(32n 2 ) is the beta function for the top Yukawa coupling, and the loop 
function 



(m£ to 2 - \ rr£ + to 2 - ( m 2 - 

depends explicitly on the renormalization scale. In the above, 6 ts = <pA t + & — <t>A t ~ <Pa s up 
to one loop accuracy determined by (30). 

3.2 The Higgs Mass Calculation 

We now turn to the Higgs mass calculation at one-loop in the presence of CP violation in 
the stop LR mixing. The mass-squared matrix of the Higgs scalars is 

Ml = (J^v) , (32) 



subject to the minimization constraints Eqs. 22-24 and (30). In the above, $j = (0j, </?;). 
Clearly, two linearly independent combinations of the pseudoscalar components (p U; d, s are the 
Goldstone bosons Gz and Gz r , which are eaten by the Z and Z' gauge bosons when they 
acquire their masses. These two modes are given by 

Gz = — sin f3ip u + cos P<pd , Gz> = cos j3 cos aip u + sin /3 cos aipd — sin af s , (33) 

and hence the orthogonal combination 

A = cos (5 sin aip u + sin (3 sin a^d + cos aip s (34) 

is the physical pseudoscalar Higgs boson in the CP-conserving limit. In the decoupling limit, 
v s ^> v , sin a — > 1 and cos a — > 0, in which case Gz and A reduce to their MSSM expressions. 
In the basis of scalars B = (fid, 0s, ^4}, the Higgs mass-squared matrix M 2 takes the form 

( M 2 U + M\ cos 2 M 2 d - M\ sin f3 cos f3 M 2 S - M\ cot a cos f3 M 2 A sin 9 t s \ 

M 2 d - M\ sin cos p Mj d + M\ sin 2 p Mj s - M\ cot a sin p Mj A sin 9 t 

M 2 S - M\ cot a cos P M 2 ds - M\ cot a sin p M 2 S + M\ cot 2 a M 2 A sin 6 t . 

\ M 2 A sin# ts M 2 A sin^ ts M 2 A sin# ts M 2 / 

in which our notation explicitly demonstrates that all of the entries M. 2 A (i — u, d, s) iden- 
tically vanish in the CP-conserving limit 6 ts — > <fio- In the above, 

Mi = M 2 (l+/^^% A (36) 
V \A S \ cosy / 



,(35) 
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which depends explicitly on the renormalization scale, and 



sin a 



+ 4/3, 



V 



m 2 |/z e// | 2 |A| 2 sin 2 6> ts 



m2 _ m 2"l sin a sin /T 



(37) 



is the one-loop pseudoscalar mass in the CP-conserving limit. The loop function Q is inde- 
pendent of the renormalization scale and has the functional form 



, \ m 2 - + m 2 - 

£?(mf ,mf) = 2--| f In 

V ti' t 2 j mA _ mA 

tl t 2 




(38) 



We now turn to the mass parameters M 2 - (i, j = u, d, s) which appear in M 2 . These entries 
may be represented as 



M. 



(47T) 



(pirn 2 + mfpj) ( 

- — r^(2 -G) + ( Pi Pj + CiO + M 



mi + 



> 4 



+ 



PiPj + 



— '2 — 2 

mfml; 



V 



[m% — m£ ) 

V tl 12/ 



G - 5 iu 5j U h t In ■ 



Q 1 



(39) 



in which \j = A^- for i ^ j, \j = 2A« for i — j. For notational purposes we have also 
introduced the dimensionless quantities 



p u = h 2 - \ u , p d = (h 2 s - \ ud )/2 , p s = (h 2 s - A us )/2 



(40) 



as well as the dimensionful ones, 



777, 



(J + h t \A t \ (\A t \ - \p eff \ cot P cos 9 t( 

,2| 



™<d = (d8 + h t \p e ff\(\p e ff\-\A t \t!mf3cos9 ts ) 



(41) 
(42) 



m c 



( s 5 + -^h 2 \p eff \ (\peff \ ~ |A t |tan/3cos0 ts ) , 



(43) 



with 5 = M~ — M~ + ( u vl + (dVd + (sV 2 - The new dimensionless couplings appearing here 
are pure D term contributions 



I/O 5 o \ 1 



— Qu)Qu9y' 



15 1 

Cd = g {gl ~ 3 0y) + g Wo ~~ Qu)Qd9Y> 

(s = -((u + (d)- 



(44) 

(45) 
(46) 



Finally, the scalar-pseudoscalar mixing entries M 2 A (i = u, d, s), which exist only if there are 
sources of CP violation in the Lagrangian (as has been made explicit in A4 2 by factoring out 
sin# ts ), are given by 



0R l/v/Pt| 

^Pht ■ o 2 

sin a mA — mA 

tl *2 



2 2 

mA — mA 

Ipi I' ? + 
mA + mA 

tl *2 



mf 



2 2 

mi — mi 

tl *2 



2 2 
Z 5 * 2 I 2~ 

mi + mA 



ti 



t2. 



(47) 
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and are scale independent. These results agree with the tree level computations of [41]. 
After identifying (11) with the |/i e //| parameter of the MSSM, the doublet sector of the 
mass-squared matrix agrees with that of the MSSM [22]. Finally, the results also agree with 
those of [45] in the CP-conserving limit (sin9 ts = 0). 

As previously stated, there are three CP even and one CP odd Higgs boson in the CP- 
conserving limit. The mass of the CP odd Higgs boson A is given in (37), while the masses 
of the CP even scalars arise from the diagonalization of the upper 3x3 subblock of (35). 
The masses and mixings then differ from their tree level values by the inclusion of radiative 
effects. In this limit, the only source of CP violation is the CKM matrix and one easily 
evades constraints from the absence of permanent EDMs for leptons and hadrons. The 
lightest Higgs boson has a larger mass than M z even at tree level, and the radiative effects 
modify it sizeably [45]. Once the radiative corrections are included a direct comparison 
with experimental results is possible. In principle, one can constrain certain portions of the 
parameter space using the post-LEP indications for a light scalar with mass > 114 GeV. 

In the presence of CP violation, there are four scalar bosons with no definite CP quantum 
number. This results from the mixing between the CP even scalars 4> u ,d,s with the CP 
odd scalar A via the entries M 2 A smd ts in (35). The main impact of the CP breaking 
Higgs mixings on the collider phenomenlogy comes via the generation of novel couplings for 
Higgs bosons which eventually modify the event rates and asymmetries. Indeed, a given 
Higgs boson can couple to both scalar and pseudoscalar fermion densities depending on the 
strength of CP violation [22]. Moreover, the coupling of the lightest Higgs to Z bosons can 
be significantly suppressed, avoiding the existing bounds from the LEP data [24, 25]. The 
CP- violating entries of Ai 2 grow with \/j, e ffA t \ as in the MSSM. The mass-squared matrix 
is diagonalized by a 4 x 4 orthonormal matrix 1Z 

n-M 2 h .n T = diag. (M 2 Hv M 2 H2 , M| 3 , M|J . (48) 

To avoid discontinuities in the eigenvalues it is convenient to adopt an ordering: Mh 1 < 
M H2 < M Hz < M Hi . The mass eigenstates H { can then be expressed as 

Hi = n iu (j) u + n id( j) d + n is( j) s + n iA A (49) 

2 m 

in which e.g. \lZiA\ is a measure of the CP odd composition of Hi. The elements of 7Z 
determine the couplings of Higgs bosons to the MSSM fermions, scalars, and gauge bosons. 

3.3 Comparison with MSSM 

Before turning to the numerical analysis, it is instructive to compare the origin of Higgs 
sector CP violation in the U(l)' models to that within the MSSM. Let us first consider the 
case of the MSSM, in which the Higgs sector consists of the two electroweak Higgs doublets 
H Utd . It is useful to start with the most general renormalizable Higgs potential for a two 
Higgs doublet model (2HDM), which must be built out of the gauge invariant combinations 
\H U \ 2 , \H d \ 2 , and H u ■ H d as follows: 

Ken DM = m 2 u \H u \ 2 + m 2 d \H d \ 2 + (mjH u -H d + h.c.) 
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+ ^i\H u \ A + A2|-£^f| 4 + A3|i/ U | 2 |-£^| 2 + \AH U ■ H d \ 2 



+ 



\ 5 (H U ■ H d f - (\ % \H d \ 2 + \ 7 \H U \ 2 ) H u -H d + h.c] , (50) 



in which m 2 , ^5,6,7 can be complex. In a general 2HDM, the Higgs sector exhibits CP 
violation if any two of these couplings have nontrivial relative phases. Spontaneous CP 
violation can also occur for certain ranges of the parameters [50] . However, at tree level the 
MSSM is a special 2HDM, with m\ = B^i = b, m 2 u d = m\ u d , and 

X 1 = A 2 = G 2 /4; A 3 = (g 2 2 - g 2 Y )/A; A 4 = -g 2 2 /2; A 5 = A 6 = A 7 = 0. (51) 

As previously discussed, there is only one complex coupling Bfj, in the MSSM Higgs potential 
at tree level, and hence its phase can always be eliminated by a suitable PQ rotation of the 
Higgs fields. Although the Higgs sector is CP-conserving at tree level, CP violation occurs 
at the loop level if Of and/or a are nonzero, with the dominant contribution involving t . If 
t 7^ 0, a relative phase between the VEVs of H u and H d is generated [22]. 

Essentially, while the U(1) P q symmetry of the MSSM forbids nonzero values of A 5)6i7 at 
tree level, these couplings are generated by radiative corrections because U(1)pq is softly 
broken by the B\x term. For example, the effective A 5 coupling which is generated at one-loop 
is approximately 

As ~ vJl ( " A ' )2; (52) 

iD7r m SUSY 

see [22] for the explicit expressions. 12 

Within the U(l)' models, the tree level Higgs potential does not allow for explicit or 
spontaneous CP violation. However, it is possible to make a stronger statement: unlike the 
MSSM, the Higgs potential in this class of U(l)' models does not allow for CP violation 
at the renormalizable level at any order in perturbation theory. To see this more clearly, 
consider the most general renormalizable Higgs potential for H u , H d , and S. The potential 
can be expressed as a function of the gauge- invariant quantities |i^«| 2 , |-£^z| 2 , \H d - H u \ 2 , and 
SH U ■ H d : 



Ken = m 2 u \H u \ 2 + m 2 d \H d \ 2 + m 2 s \S\ 2 + (m 12 SH u ■ H d + h.c.) 
+ ^u\H u \ 4 + \ d \H d \ 4 + \ S \S\ A 

+ ^ud\H u \ 2 \H d \ 2 + \ US \H U \ 2 \S\ 2 + \ ds \H d \ 2 \S\ 2 + \ ud \H u ■ H d \ 2 , (53) 

At tree level, the dimensionful parameters m 2 d = m 2 Hud and myi = h s A s , and the dimen- 

sionless couplings have all been listed before except \ ud = \g\ — h 2 s . Therefore, even in the 
most general renormalizable Higgs potential there is only one coupling which can be complex 
(777.12); this is because the gauge-invariant operator SH U • H d is already dimension 3. Hence, 
the global phases of the Higgs fields (more precisely of the combination SH u H d ) can always 

12 Note that spontaneous CP violation (SCPV) requires m§ < \^v u Vd- As A5 is loop suppressed in the 
MSSM, SCPV would require a very small m§, leading to an unacceptably light pseudoscalar Higgs mass [50]. 
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be chosen such that the phase of mi 2 is absorbed. Note that this statement, while true for 
the tree- level potential of (6), does not depend in any way on perturbation theory. 13 

As the Higgs potential conserves CP to all orders at the renormalizable level, CP violation 
can enter the theory only through loop-induced nonrenormalizable operators. The form of 
(17) demonstrates that the one-loop contributions to the Higgs potential include an infinite 
series of terms involving powers of the Higgs fields. While these terms include contributions 
to the potential at the renormalizable level, they also include a tower of nonrenormalizable 
terms, such as 

V nr = ...+ (^^(SH U ■ H d f + h.c) + . . . , (54) 

\ m SUSY / 

in which m SUSY denotes a typical sfermion mass. By U(1)r invariance, the coupling A of 
the (SH U ■ H d ) 2 term is proportional to A ~ A 2 / (lQir 2 m 2 SUSY ) ■ Such a term is generated by 
the one-loop diagram formed from the Lagrangian interactions h s h*SH^u* L u c ^ + h.c. (from 
F terms) and the soft SUSY breaking interaction h t A t uiju c L + h.c. For (S) ^> (H Utd ), (54) 
effectively leads to the coupling 14 

-$^(H U -H d ) 2 , (55) 

m SUSY 

with 



ie// ^ (A t (s) ) 2 

lQn 2 m 2 



w - • (56) 

l SUSY 

In general, one can expand the one-loop potential in powers of the phase-sensitive gauge- 
invariant operator SH U • H d : 



5V = ... - (3 ht h s T (Q 2 , m|, m|) A t SH U ■ H d 



+ (3 ht h 2 h 2 s Q (mf^mfj (SH U ■ H d f + h.c. + . . . , (57) 

( — m£ ) 

\ tl t-2/ 

in which we have presented only the phase-sensitive corrections up to quadratic order (this 
expansion can of course be continued to higher orders with no difficulty at all). The first term 
renormalizes the h s A s SH u ■ H d operator in the tree level potential, while the second term is a 
new higher-dimensional operator. Both terms violate CP through the phase of A t (S) (recall 
this phase is irremovable if A t and A s have a nontrivial relative phase 9 ts ). The effective 
theory at scales below (S) is equivalent to the MSSM (with /i and Bfx parameters related to 
the other soft parameters of the model). One concludes from (57) that the size of the CP 
violation in the Higgs sector depends on the extent to which the U(l)' breaking scale is split 
from the electroweak scale. Indeed, below the scale (S) , the coefficients of the CP- violating 
effective operators in (50) grow with \A t \v s (or equivalently \A t \ \/i e ff\), in agreement with 
the CP- violating M 2 U d ^ A entries of the Higgs mass-squared matrix. 

13 Note that the structure of the potential is very different in the case of the NMSSM, in which the S is 
a total gauge singlet. As gauge invariance then does not restrict the possible S couplings, the Higgs sector 
generically violates CP at tree level [51]. 

14 Note that SCPV is also not viable in this potential, for the same reason as in the MSSM. 
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4 Phenomenological Implications 



In this section we discuss the existing constraints on U(l)' models as well as their phe- 
nomenological implications with explicit CP violation. 

4.1 Constraints from Z — Z' Mixing 

In the previous section, we computed the radiatively corrected Higgs boson mass-squared 
matrix (35). If the eigenvalues of the Higgs mass-squared are all positive definite, the pa- 
rameter space under concern corresponds to a minimum of the potential. The parameter 
space is of course also constrained by the fact that direct collider searches have yielded lower 
bounds on the sparticle and Higgs masses. Within U(l)' models, further constraints arise 
from the nonobservation to date of a Z', both from direct searches [38] and indirect precision 
tests from Z pole, LEP II and neutral weak current data [28, 39]. The strongest constraints 
arise from the mixing mass term between the Z and the Z' induced by electroweak breaking: 

Mf A 2 \ 

(58) 



M z -z' 
in which 



Mi, J 



Mf = G/V/4 (59) 
Ml = 9 2 Y/ (Qy u + Q 2 d v 2 d + Q 2 s v!) (60) 

A 2 = l -g Y ,G(Q u vl-Q d v 2 d ). (61) 

Current data requires A 2 <C Mf,, Mf, because the Z — Z' mixing angle 

1 / 2A 2 \ 

^' = 2 arCtan U 2 ,-M 2 J- (62) 

must not exceed a few x 1CT 3 in typical models. 

Let us review the implications of this constraint, which was studied in [28, 41]. One 
can see from (62) that unless M Z ' 3> M z , the Z — Z' mixing angle is naturally of 0(1). 
Therefore, a small az-z' requires a cancellation in the mixing term A 2 for a given value 
of tan/5. For models in which Mz> ~ O(Mz), this cancellation must be nearly exact; this 
can be slightly alleviated if the Z' mass is near its natural upper limit of a few TeV. Hence, 
tan 2 /5 must be tuned around Qd/Q u with a precision determined by the size of az-z'- In 
our analysis, we will eliminate tan/3 from (62) for a given value of az~z r - 



tan2 , vQt-oz-z'i-l+rfiQS + Q, 2 ^ 

vQu + az-z' (-i+v 2 {Qu 2 + Q s 2 $y 



in which r] = 2gy>/G, and we used tan(2 otz-z 1 ) ~ ^olz-z 1 - Having fixed tan/5 in this way, 
a multitude of parameters remain which can be varied continuusly as long as all collider 
constraints are satisfied. In [41], two phenomenologically viable scenarios were identified: 
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• Light Z' Scenario. Clearly, the U(l)' symmetry can be broken along with the SM 



gauge symmetries at the electroweak scale. 15 In this case v s ~ v, tan/5 ~ J\Qd\/\Q- 



and Mz> is of order Mz (the precise factor depends on the size of gY>\Q s \). However, 
the collider constraints on such a light Z' are severe within typical models, and hence 
it can be accommodated in the spectrum only if it is sufficiently leptophobic. Note 
that within this framework, leptophobic U(\)' couplings lead to a generic difficulty 
related to lepton mass generation: as Qu A 7^ 0, if the electron mass is induced via the 
Yukawa coupling h e LiHdEf, the leptons necessarily have nonvanishing U(l)' charges. 
The electron mass (and perhaps all light fermion masses) then must be generated via 
nonrenormalizable interactions which guarantee the neutrality of L\ and E\ under the 
£7(1)'. In practice, this would need to be investigated within specific models. 16 

• Heavy Z' Scenario. In this scenario, the U(l)' breaking is radiative (driven by the 
running of m| to negative values in the infrared) and occurs at a hierarchically larger 
scale than the electroweak scale. However, gauge invariance does not allow for the U(l)' 
and electroweak breakings to decouple completely (as Qu ud 7^ 0). The electroweak 
scale is then achieved by a cancellation among the soft masses, which are typically of 
O(Mz'), with a fine-tuning 0(Mz> / Mz). As discussed in [41], excessive fine tuning 
is avoided if Mz> in units of the heavy scale is roughly bounded by the ratio of the 
charges, Min[|<5 s /<9d|, \Qs/Qd\\- There are several advantages of the heavy Z' scenario. 
First, the Z — Z' mixing can be kept small enough with less fine-tuning of the A 2 
in (58); in particular, Q u = Q,i is no longer a requirement. In addition, the collider 
constraints are less severe for Z' bosons with TeV-scale masses in typical models; for 
example, leptophobic couplings are not generically a phenomenological necessity. 

4.2 Constraints from Dipole Moments 

Let us now turn to dipole moment constraints. Recall in SUSY theories dipole moments 
of the fundamental fermions are generated by gaugino/Higgsino exchanges accompanied by 
sfermions of the appropriate flavor. The dipole moment under concern may (e.g. the electric 
and chromoelectric dipole moments of the quarks) or may not (e.g. the anomalous magnetic 
moment of muon) require explicit sources of CP violation. 

15 As shown in [41], at tree level a light Z' boson with a vanishing Z — Z' mixing (for Q u — Qd) naturally 
arises when \A S \ is the dominant soft mass in the Higgs potential. Such trilinear coupling induced minima can 
also accommodate a heavy Z' boson. This can happen in models in which there are several additional singlets 
in a secluded sector coupled to the Higgs fields H u ^ and S via the gauge or gravitational interactions [44]. 
Furthermore, these large trilinear coupling scenarios (with light Z' bosons) also have interesting implications 
for baryogenesis, due to the first order phase transition at tree level. If the phase transition remains first 
order after radiative corrections are included, then 9 may be sufficient to generate the baryon asymmetry. 
The electroweak phase transition in Z' models with a secluded sector is strongly first order (with a heavy 
enough Z' without any fine-tuning), and electroweak baryogenesis in such models can be viable in a greater 
region of parameter space than in the minimal model [52]. 

16 However, the kinetic mixing between the hypercharge and Z' gauge bosons can be used to decouple 
leptons from Z' though all leptons, with nonzero U(l)' charges, acquire their masses from their Yukawa 
couplings [48]. 
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In the MSSM, dipole moments can provide important constraints on the parameter space. 
For example, the anomalous magnetic moment of the muon is in principle an important 
observable either for discovering SUSY indirectly or constraining SUSY parameter space; 
however, at present the theoretical uncertainties present in certain nonperturbative SM con- 
tributions lead to difficulties in carrying out this procedure using recent data (see e.g. [53] for 
a review of the basic physics and [54] for the most recent experimental results). At present, 
the most stringent constraints arise from electric dipole moments (EDMs). As is well known, 
the experimental upper bounds on the EDMs of the electron, neutron, and certain atoms 
impose particularly severe constraints on the parameter space of general SUSY models. In 
contrast to the SM, in which EDMs are generated only at three-loop order (as the only 
source of CP violation is in the CKM matrix), the sources of explicit CP violation in SUSY 
theories include phases in flavor-conserving couplings which, if present, lead to nonvanishing 
one-loop contributions to the EDMs which can exceed the experimental bounds. As these 
phases generically filter into the Higgs sector, it is important to include the parameter space 
constraints provided by the EDM bounds. 

Let us consider the dipole moments which arise within this class of U(l)' models. After 
replacing the \x parameter of the MSSM by fi e ff in (11), all one-loop dipole moments are 
found to be identical to their MSSM counterparts except for an additional contribution 
generated by the Z'-f diagram in Fig.l. Here Z' is the U(l)' gaugino with mass My. This 
diagram generates the operator Dff L a^ u F^f R , in which 

D,M ~ 2^ [1-4/1 *«♦.'.-*.> - R,\M.f\ . (64) 

/ 

in units of the electromagnetic or strong coupling. In the above Rf = (tan (3)~ 2If , My char- 
acterizes the typical sfermion mass, and recall that the reparameterization invariant phases 
are defined in (4). Clearly, the (chromo-)electric and (chro mo-) magnetic dipole moments are 
generated, respectively, by Im[Dj] and Re[.D/]. The expression above is approximate esti- 
mate (valid in the limit that Mj~^> My) of the exact amplitude; a more precise treatment 
would take into account the mixing of all six neutral fermions. The amplitude (64) is similar 
to the bino exchange contribution in the MSSM. 

Within the aforementioned light Z' scenario, for phenomenologically viable models D e (Z') 
vanishes because the lepton couplings to the Z' are necessarily leptophobic. Therefore, for 
instance, the electron EDM is completely decoupled from the presence of an electroweak 
scale U(l)' symmetry. This conclusion extends to other leptons for family universal Z 1 
models. This may also be relevant for the hadronic dipole moments depending on whether 
or not the Z 1 boson is hadrophobic (assuming it is detected in present and/or forthcoming 
colliders). As |/i e //| <^ Mj within the light Z' scenario, the dipole moments of both the 
up-type and down-type fermions are largely controlled by the corresponding Af parameters. 
In contrast, the U(l)' charges are not necessarily suppressed for any fermion flavor in the 
heavy Z' scenario and thus the Df(Z') contribution to dipole moments can compete with the 
MSSM amplitudes. In this scenario, |// e //| ~ Mz> Mz, and hence both terms in Df(Z') 
are important. The dipole moments become sensitive to 9y s + 9 if the Af parameters are 
sufficiently small compared to \fJ- e ff\- 
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Figure 1: The Z'-sfermion diagram which contributes to the (chromo-)electric and (chromo-)magnetic 
dipole moments of the fermion /. The photon (7) or gluon (g) are to be attached in all possible ways. 

As the dipole moments genetically scale as mf/M~, (which is clear from the form of 
(64)), when Mj~ 0(Mz) the EDMs typically exceed the existing bounds by 2 to 3 orders of 
magnitude if the phases are 0(1). As discussed briefly in the Introduction, one possibility is 
satisfying the experimental bounds while retaining 0(1) phases is to raise Mjto multi-TeV 
values, which in effect requires the sfermions of first and second generations to be ultraheavy 
[2, 4]. Another way of suppressing the EDMs is to invoke accidental cancellations between 
different contributions, i. e. to find regions of parameter space in which the SUSY amplitudes 
interfere destructively. In the MSSM with low values of tan /3 , this has been shown to occur 
with almost no costraint on any of the invariant phases except \9~±\ — |</y + 4>m 2 \ ^ 7r/10 
[4, 5, 6, 7, 8]. This strong constraint follows from the fact that gy <C #2, and thus the SU(2) 
gauginos dominate the EDMs. Within the U(l)' framework, the EDM constraints can have 
varying implications depending on the size of (64). 

• If gy ~ 0(gy) or (more generally) gy <C #2, the EDM constraints on the parameter 
space are similar to that of the MSSM except for a slight folding of the cancellation 
domain due to the inclusion of (64). Once again, the most strongly constrained phase 
is 9~±, with \9~±\ < 7r/10 in the low tan/5 regime. As 9~± = 02 S + 9 and 9 is a 
loop-suppressed angle (30), the EDMs provide a constraint on 9 2s '- \02s\ ~ tt/10. Con- 
sequently, the dynamical solution to the \i problem present in this class of U(l)' models 
also solves the SUSY CP hierarchy problem in specific models of the soft parameters 
in which (at least) the SU (2) gaugino mass has the same phase as the A s parameter 
(then 9 2s = by definition). 17 

• If gy> ^ g-ii the dipole moment amplitude Df(Z') becomes comparable to or larger than 

17 In this paper, we have not addressed the origin of the phases of the soft parameters in (3), and hence 
we cannot make any claims about how one solves the SUSY CP hierarchy problem within this framework. 
However, it is worthwhile to note that models of the soft parameters in which the gaugino masses and A 
terms have the same phases are quite common within various classes of four-dimensional string models (at 
least at tree level) under plausible assumptions [55]. 
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the SU (2) gaugino contribution, and the cancellation domain found in the MSSM will 
be significantly folded. In this case, the EDMs will constrain a combination of the 
phases in (64) and 9~±. Such a scenario, however, can have tension with the standard 
picture of gauge coupling unification at a high fundamental scale (although in principle 
it could be considered as a possibility in generic low scale realizations). 

Until this point, we have only discussed one- loop EDMs. It was pointed out a while ago 
[11] that in certain regions of MSSM parameter space certain two-loop contribution contri- 
butions which exclusively depend on the third generation sfermions can be nonnegligible. 
These contributions, which are particularly relevant if the one-loop EDMs are suppressed 
solely by ultraheavy first and second generation sfermion masses, involve the same phases 
which predominantly filter into the Higgs potential at one-loop (i.e. the phases present in 
the stop mass-squared matrix). However, these two- loop EDMs become sizeable only at 
large tan/3. In this paper, we have restricted our attention to small tan/5 values, which 
is a well- motivated parameter regime (e.g., tan/3 = 1 is allowed within this framework, in 
contrast to the MSSM). Hence, these contributions will not provide significant parameter 
space constraints in our numerical analysis. 

4.3 Numerical Estimates for Higgs Sector CP Violation 

In this section, we present sample numerical calculations of the Higgs boson masses and 
mixings derived in Section 3.2, taking into account the phenomenological constraints on the 
parameter space discussed in Sections 4.1 and 4.2. 

In the absence of CP violation, the scalar-pseudoscalar mixing terms of the Higgs mass- 
squared matrix (35) vanish (sin (9 = 0), and hence (35) takes on a block diagonal form. The 
structure of (35) demonstrates that in this limit there is one CP even scalar with mass oc v s 



and a CP odd scalar with mass proportional to y \A s \v s . In addition, there is a light CP even 
scalar of mass ~ Mz and a heavier CP even scalar with its mass controlled by a combination 
of v and M A . However, in the presence of explicit CP violation, the Higgs bosons cease to 
have definite CP parities. The strength of CP violation in the Higgs sector is parameterized 
by the reparameterization invariant phase 6 ts , which induces a nonvanishing through the 
relation (30). The induced phase 9 is a loop- induced and scale-dependent quantity which is 
particularly enhanced in parameter regions with a low Ma- 

As discussed in Section 4.2, while the one-loop EDM constraints strongly constrain the 
phase 02 S , this phase is not the dominant source of CP violation in the Higgs sector for 
small values of tan (3 and hence this constraint does not restrict the parameter space for our 
analysis. The dominant corrections to the Higgs potential arise from top and stop loops, 
and the dipole moments of the fermions in first two generations feel such effects only at two 
loop level. In fact, in low tan f3 limit (which is the domain in which our analysis of the Higgs 
potential is valid), such effects are completely negligible [11]. Therefore, the EDM constraints 
do not have a direct impact on our analysis of CP violation in the Higgs sector (we simply 
assume that the dipole moment constraints have been saturated either via cancellations or 
by choosing the first and second generation sfermion masses heavy enough; we could also 
simply assume that all phases except 9 ts are small). 
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We now turn to the analysis of the parameter space, including the nontrivial constraints 
arising from Z — Z' mixing. The fundamental parameters relevant for the Higgs sector include 
{v s , A s , A t , Mq, My,h s , Q u , Qd, gy, &ts\- We fix a subset of these parameters as follows: (i) 
olz-z 1 = 10~ 3 , which is well below the present bounds; (ii) g\, = (5/3)G 2 sin 2 9w, as inspired 
from one-step GUT breaking; (iii) h s = l/y/2, as motivated by the RGE analysis of [41]; (iv) 
Qu = Qd = — 1, such that tan/3 remains close to unity (as can be seen from (63)); and finally 
(v) Mq = My. The remaining parameters can be fixed on a case by case basis depending 
on the range of values assumed for Mz>- A few notational comments are also in order. 
Although (35) suggests that M\ can be chosen to be a fundamental parameter and this is 
what is traditionally done in the MSSM, we prefer to work instead with A s for consistency 
with previous discussions in this paper as well as the tree level analysis of [41]. In addition, 
in our numerical results we fix the renormalization scale to be Q = (2m t + M z >)/2. This 
differs once again from the MSSM, where the renormalization scale is chosen to be Q = m t 
in order to minimize the next-to-leading order corrections. Such higher order corrections are 
beyond the scope of this paper; our choice for Q can be regarded as some nominal value in 
between the electroweak and U(l)' breaking scales. 

We begin with an analysis of the light Z' scenario. For purposes of definiteness, we set 
Mq = 2v s , v s = v/y/2 ~ m t , and \A S \ = v s , in which case Mz> — 2M Z and \(J> e ff \ — M z . The 
SUSY phase 9 ts influences both the Higgs masses and their mixings, as shown in Figure 2. In 
the left panel, the variation of the lightest Higgs mass with 9 ts is displayed for several values 
of \A t \/v s . For \A t \/v s = 1/2, 1 and 2 M Hl grows gradually with 9 ts , peaking at 9 ts = it. 
This behaviour is easy to understand: as the magnitude of the stop LR mixing depends 
strongly on 9 ts , the variation of Mh 1 with respect to 9 ts simply displays the well-known fact 
that the lightest Higgs mass depends strongly on the value of the stop mixing. Indeed, 

\ M LR\e t s=K = \M + |/A e// |cot/3 
\ M LR\e ts =i) \A t \ - |/v/|cot/3' 

which becomes large when \A t \ and \fJ> e ff\ are of comparable size. The ratio (65) gets sat- 
urated with further increase of \A t \; however, in this case |A t ||// e ^| also becomes large, 
which affects both the Mp and Mf A entries of the Higgs mass-squared matrix. While the 
former shifts the peak value of Mh 1 towards the point of maximal CP violation (see the 
dot-dashed curve in the figure), the latter enhances the scalar-pseudoscalar mixings. The 
generic strength of the scalar-pseudoscalar mixings can be determined e.g. by working out 
the CP-odd composition of H 3 (the would-be pseudoscalar Higgs). The result is shown in the 
right panel of Figure 2. Clearly, the M? A sm0 ts elements of the Higgs mass-squared matrix 
are not large enough to enhance such mixings (|7?.3a| 2 falls at most to 99.75% for \A t \ = 4i> s ). 

The functional dependence of the heavier Higgs boson masses on 9 ts is opposite that 
of M Hl in that the masses tend to decrease as 9 ts ranges from to n; e.g. when \A t \ = 
Av s , (M H4 , M H3 , M H2 ) fall from (245,224,191) to (234,206,182) GeV. In accord with the 
analytical expression (30), 9 grows with \A t \ until it arrives at the peak value of ~ 30% for 
\A t \ — Av s for maximal CP violation. For low Mz> minima, the scalar-pseudoscalar mixings 
(which govern the novel CP violating effects in the Higgs couplings to fermions, gauge bosons 
and other Higgs bosons) are typically small due to the low value of |/i e //|- 



20 



> 

CD 



136 
134 



132 - 



130 - 



128 



i — i — i — i — i — r 




j i i i i L 



0.5 1 1.5 2 2.5 3 

&ts 



100 



99.95 - 
99.9 - 



« 99.85 - \ 



99.8 h 
99.75 



\\ 

V 



s ■ , 



/ / 
/ / 
/ / 



j i y i i l 



0.5 1 1.5 2 2.5 3 

&ts 



Figure 2: The 6 ts dependence of the lightest Higgs mass and the CP-odd composition of H 3 in the light Z' 
scenario. The solid, dashed, dotted, and dot-dashed curves correspond to, respectively, |A t |/i> s = 1/2,1,2 
and 4 with v s = v/\/2. 




Figure 3: The ts dependence of the lightest Higgs mass (left panel) and the CP-odd composition of H3 
in the heavy Z' scenario. Here solid, dashed, dotted, and dot-dashed curves correspond, respectively, to 
\A s \/v s = 1/5,1/2,3/4 and 1 with v s = 1 TcV. 
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We now discuss the heavy Z' scenario, setting v s — 1 TeV, M~ = 750 GeV, and \A t \ — 
2Mq. Figure 3 depicts the variations of the lightest Higgs mass (left panel) and the CP-odd 
composition of the would-be Higgs scalar as a function of 9 ts and \A S \. In both figures, 
the solid, dashed, dotted, and dot-dashed curves correspond, respectively, to |Ag|/i> s = 
1/5,1/2,3/4 and 1. In contrast to the light Z' scenario as shown in Figure 2, here we 
illustrate the dependence on \A S \ (or equivalently Ma), as this parameter remains largely 
free in the heavy Z' limit [41]. As the left panel of the figure shows, the mass of the lightest 
Higgs is typically larger than that in the light Z' scenario. The lightest Higgs mass is also a 
steep function of \A S \, which becomes increasingly smaller as Ma increases due to decoupling. 
Note also that the dependence of M Hl on 6 ts in this scenario is similar to the case within the 
light Z' scenario; once again, this is because the radiative corrections to the lightest Higgs 
mass strongly depend on the value of the stop mixing parameter. 

However, in contrast to the light Z' scenario, the scalar-pseudoscalar mixings in the 
heavy Mz> limit are sizeable, as shown in the right panel of Figure 3. This feature is ex- 
pected because the strength of the Higgs sector CP violation is governed by the size of the 
singlet VEV, i.e. the effective /i parameter, and in this scenario |/i e //| ~ Mz>- In general, 
the CP- violating mixings grow larger as A s decreases, because in this case Mf A sin 9 t s can be 
comparable to Ma, which facilitates scalar-pseudoscalar transitions. For \A S \ = v s /5, the CP- 
odd composition of the would-be pseudoscalar Higgs falls down to 70% around $t s ~ 7r/6. 
However, as \A S \ increases (while keeping |/i e //| and \A t \ fixed), the diagonal elements of 
(35) also increase, with the result that the CP-violating effects become weaker. The large 
variations in 17^3^ | 2 depicted here are due to the mixings between H 3 and H 2 . Indeed, for 
ILIA'S — 1/5, 1/2, 3/4 and 1 the two masses are strongly degenerate, with (M# 2 , Mh 3 ) start- 
ing at (476,477), (722,726), (876,881), (1013, 1007) and decreasing to (417,418), (685,688), 
(846,851), (987,981) GeV as 9 ts varies from to n. Note that the scalar-pseudoscalar 
conversions are more efficient when the two masses are highly degenerate. 

For the values of | A s | / v s exhibited above, the Higgs sector is within the decoupling regime 
(Ma > 2Mz) 18 , in which the lightest Higgs resembles the SM Higgs boson, the heaviest Higgs 
is singlet-dominated with a mass of order Mz', and the two intermediate mass Higgs (the 
CP odd scalar and the second heaviest CP even scalar in the absence of CP violation) 
are strongly degenerate. The lightest Higgs boson is essentially CP even (l-Ri^l 2 <C 0.1% for 
IAEA's > 1/15) and hence is decoupled from CP- violating effects, although its mass depends 
strongly on Q ts . However, there are phenomeno logically interesting corners of parameter 
space with sufficiently small values of |A,|/t> s in which the lightest Higgs boson can have 
a significant mixing with the would-be pseudoscalar. As the lightest Higgs mass is a steep 
function of \A s \/v s , for a value of M Hl consistent with LEP bounds the CP-odd composition 
of Hi cannot be larger than 20%. It is important to keep in mind though that the couplings 
of the lightest Higgs boson to gauge bosons and fermions are modified when the lightest 
Higgs has a significant mixing with the would-be pseudoscalar (the modifications grow with 
the CP-odd composition of the lightest Higgs), such that the existing LEP bounds may not 
be applicable (see e.g. [24, 25] for discussions within the MSSM). 

18 See [22, 23] for a more precise definition of the decoupling regime in the CP-violating MSSM. 
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Our results demonstrate that the CP-violating effects in the Higgs sector, or more pre- 
cisely, the mixing between the would-be scalars and pseudoscalars in the CP conserving 
limit, are generically highly suppressed in the light Z' models but can be sizeable in the 
heavy Z' scenario, even though the masses can vary strongly with 6 ts (which is of course a 
CP-conserving effect). This behaviour is exactly in accordance with the general discussion 
of Section 3.3, in which we demonstrated that the CP-violating terms in the Higgs potential 
necessarily originate from nonrenormalizable terms present at one-loop (such terms are en- 
coded within the full Coleman- Weinberg potential). The strength of such terms in e.g. the 
doublet sector then scale according to the ratio of the singlet VEV v s ~ \/J> e ff\ to the scale of 
a typical soft mass. Hence, within the light Z' scenario (in which the effective \i parameter 
is small) CP- violating effects are suppressed, while the large |/i e //| present in the heavy Z' 
scenario can allow for spectacular effects of CP violation. 

We close this section with a brief discussion of the implications for collider searches. In 
general, at least a subset of the Higgs masses within this class of U(l)' models can be observ- 
able at forthcoming colliders and future colliders such as TESLA and NLC. next generation 
of colliders such as TESLA and NLC. Within light Z' models, all of the Higgs bosons remain 
light after including radiative corrections, but such models generically have very small CP- 
violating Higgs couplings. In contrast, large Z' models can have large CP- violating Higgs 
couplings. As the viable regions of space typically correspond to the decoupling limit in 
which all of the Higgs bosons except the lightest Higgs are heavy, detecting the CP- violating 
effects within this scenario is similar to that within the MSSM for a large /i parameter. Such 
effects have been studied in [22, 24, 25], where it is known that Higgs sector CP violation 
can introduce sizeable modifications in the couplings of the Higgs bosons to fermions and 
vector bosons, and strongly affect the bounds inferred from the CP-conserving theory. Fur- 
thermore, the CP purity of the Higgs bosons (assuming that the collider searches establish 
their existence) can be tested by measuring CP violation in its decays into heavy quarks or 
vector bosons [56]. 

5 Summary 

In this paper, we have discussed the nature and implications of explicit CP violating phases 
present in the soft breaking Lagrangian within supersymmetric models with an additional 
U(l) gauge symmetry and an additional SM gauge singlet S. This class of models is worthy of 
further study not only because such gauge extensions are ubiquitous within four- dimensional 
string models (and other plausible extensions of the MSSM), but also they provide an elegant 
framework in which the fi problem of the MSSM. The solution, which is to forbid the bare \i 
term by U(l)' gauge invariance and generate an effective fi parameter through the VEV of 
the singlet S, is similar to that found within the NMSSM (but its generic cosmological and 
CP problems). Our results can be summarized as follows: 

• All reparameterization invariant phases can be expressed as linear combinations of 
0f s = (pA f — (ft a s and 9 as = (pM a — <Pa s (and hence a "natural" basis can be obtained 
by using U(1) R to set <p As = 0). 
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• The Higgs sector is manifestly CP conserving at tree level (and indeed at renormalizable 
level to all orders in perturbation theory). However, the CP- violating phases present 
in the stop mass-squared matrix filter into the Higgs sector at the nonrenormizaliz- 
able level at one-loop. The CP-violating effects are particularly enhanced when U(l)' 
symmetry is broken near the sparticle thresholds. 

• The spontaneous breakdown of the U(l)' symmetry near the weak scale stabilizes not 
only the modulus of fx but also its phase. The phase of fx itself is of course not a 
basis-independent quantity; however, in the "natural" basis defined above, this phase 
{9 in this basis) arises only at the loop level and is typically 1-10%, depending on the 
size of Ma (the pseudoscalar Higgs boson mass in the CP conserving limit). 

• The absence of permanent EDMs for leptons and hadrons (even assuming either can- 
cellations and/or heavy first and second generation sfermions) strongly bound the 
reparameterization invariant phase present in the chargino mass matrix (0^ + <Pm 2 ) = 
(9 — (pA s + 0m 2 )j while the other SUSY phases remain largely unconstrained. In spe- 
cific models in which the phase difference between (at least the SU(2)) gaugino mass 
parameters and A s is vanishingly small, this "SUSY CP hierarchy problem" is resolved 
because the radiative phase 9 is sufficiently small to be easily allowed by EDM bounds. 

• The CP-violating effects in the Higgs sector are quite distinct for the two phenomeno- 
logically viable scenarios with acceptably small Z — Z' mixing, because these effects 
are proprotional to the size of the effective fx term. In scenarios with a light Z 7 , 
CP-violating effects are suppressed, while heavy Z 1 models can exhibit significant CP- 
violating scalar-pseudoscalar mixings, with phenomenological implications similar to 
that of the MSSM with large fx parameter. 
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